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ABSTRACT 

We compute the cross-correlation between the Warm-Hot Intergalactic Medium and maps of cosmic 
microwave background temperature anisotropics using a log-normal probability density function to 
describe the weakly non-linear matter density field. We search for this contribution in the data mea- 
sured by the Wilkinson Microwave Anisotropy Probe. We use a template of projected matter density 
reconstructed from the Two-Micron All-Sky Redshift Survey as a tracer of the electron distribution. 
The spatial distribution of filaments is modeled using the recently developed Augmented Lagrangian 
Perturbation Theory. On the scales considered here, the reconstructed density field is very well de- 
scribed by the assumed log-normal distribution function. We predict that the cross-correlation will 
have an amplitude of 0.03 — 0.3/iK. The measured value is close to 1.5/iK, compatible with random 
alignments between structure in the template and in the temperature anisotropy data. Using the Wl 
Differencing Assembly to remove this systematic gives a residual correlation dominated by Galactic 
foregrounds. Planck could detect the Warm-Hot Medium if it is well traced by the density field re- 
constructed from galaxy surveys. The 217GHz channel will allow to eliminate spurious contributions 
and its large frequency coverage can show the sign change from the Rayleigh- Jeans to the Wien part 
of the spectrum, characteristic of the thermal Sunyaev-Zeldovich effect. 

Subject headings: cosmic background radiation - cosmology: observations - cosmology: theory - inter- 
galactic medium 



1. INTRODUCTION. 

The baryon budget in the local universe shows a deficit 
relative to the predicted density synthesized in the big- 
bang (Fukugita et al. 1998; Fukugita & Peebles, 2004). 
Galaxies and clusters contain about 10% of the total 
number of baryons and an extra 5% could be in the 
form of circumgalactic medium (CGM) around galax- 
ies, although the results of Gupta et al. (2012) of a 
large-scale massive hot gaseous halo around the Galaxy 
have been disputed by Wang & Yao (2012). Of the 
reminder 85-90%, only half has been accounted for in 
the low redshift intergalactic medium (Danforth & Shull 
2008). Hydrodynamical simulations predict that the rest 
could reside within mildly-nonlinear structures with tem- 
peratures 0.01 — IKeV, called Warm-Hot Intergalactic 
Medium (WHIM). The baryon fraction in this medium 
could be 40% (Cen & Ostriker 1999; Dave et al. 1999, 
2001, Smith et al. 2011). In the X-ray, the WHIM sig- 
nature has been searched both in emission and in ab- 
sorption. Soltan (2006) looked for the extended soft X- 
ray emission around field galaxies but his task was com- 
plicated by the need to subtract all systematic effects 
that could mimic the diffuse signal. The recent obser- 
vational effort has concentrated in searching for absorp- 
tion lines due to highly ionized heavy elements from the 
far-ultraviolet to the soft X-ray (see Shull et al. 2012, 
for a review). Alternatively, as the WHIM is highly 
ionized, in Atrio-Barandela & Mucket (2006) and Atrio- 
Barandela et al. (2008) we suggested that it would gen- 
erate measurable temperature anisotropics on the cos- 
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mic microwave background (CMB) due to the thermal 
and kinematic Sunyaev-Zeldovich effect (TSZ and KSZ, 
respectively; Sunyaev & Zeldovich 1970, Sunyaev & Zel- 
dovich 1972). Our expectations were confirmed by Hall- 
man et al. (2007) who found, using numerical simula- 
tions, that after the contribution of resolved clusters is 
removed, about one-third of the SZ flux from unresolved 
sources would be generated by unbound gas. More re- 
cently, Lieu & Duan (2013) suggested that the line-of- 
sight column density of the ionized baryons in the local 
universe could be determined by monitoring quasar light 
curves. 

Our model assumes that the undetected baryon phase, 
in the form of filaments of hot and low density IGM, 
is well described by a log-normal probability distribu- 
tion function. The filamentary structure of the inter- 
cluster medium has been recently confirmed using ob- 
servations on interacting clusters. A joint analysis of 
ROSAT X-ray and Planck CMB data has provided the 
first detection of hot and diffuse intercluster gas, extend- 
ing beyond the virial region of the cluster pair A399- 
A401 (Planck Collaboration, 2013). If the WHIM can 
be traced with templates constructed from the matter 
distribution, then cross-correlating those templates with 
CMB maps could provide further evidence of the ex- 
istence of the WHIM or provide strong constraints on 
its spatial distribution. Based on our log-normal model, 
in this article we compute the cross-correlation between 
CMB data with a template that traces the WHIM dis- 
tribution. If galaxies from the Two-Micron All-Sky Red- 
shift Survey (2MRS) are good tracers of intergalactic 
gas, we can predict theoretically the amplitude of the 
TSZ component due to the WHIM. We compare our pre- 
diction with the cross-correlation derived from Wilkin- 
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son Microwave Anisotropy Probe (WMAP) data (Jarosik 
et al. 2011). We show that the cross-correlation of 
galaxy templates and CMB data is dominated by ran- 
dom alignments. The differences between channels are 
due to Galactic foreground residuals and instrumental 
noise. Planck could provide the first detection of the 
WHIM with our method by using the 217GHz channel 
to remove systematics. Briefly, in Section 2 we describe 
our model; in Section 3 we describe the construction of 
the matter templates, show that they are well describe by 
our log-normal formalism and compute the theoretically 
expected amplitude of the cross-correlation; in Section 4 
we particularize our methods to a matter reconstruction 
of the 2MRS catalog and WMAP data and discuss the 
prospects of a successful detection using the forthcoming 
Planck data; finally, in Section 5 we present our main 
conclusions. 

2. A LOG-NORMAL MODEL OF THE WARM HOT 
INTERGAL ACTIO MEDIUM. 

In our model, baryons in the WHIM are distributed like 
a log-normal random field. The log-normal distribution 
was introduced by Coles & Jones (1991) to describe the 
non-linear distribution of matter in the universe. In this 
approximation, the number density of baryons ns(x, z) 
is given by (Choudhury et al. 2001, Atrio-Barandela & 
Miicket 2006) 



n B (x, z) = n (z)e 



Mx,z)-A|( Z )/2 



(1) 



where x denotes the spatial position at redshift z and 
|x(z)| is the proper distance, A^(z) =< 5 B (x, z) >, with 
<5b(x, z) is the baryon (linear) density contrast, n$(z) = 
Pb(1+z) 3 I iiBm p with ps, m p are the baryon density and 
the proton mass, respectively; finally, fis = 4/(8 + 5Y) is 
the mean molecular weight of the IGM and Y = 0.24 is 
the helium weight fraction. The log-normal distribution 
has been found to describe very well the matter statistics 
at scales larger than 7/i~ 1 Mpc based on the improved 
Wiener density reconstruction from the Sloan Digital Sky 
Survey (see Kitaura et al. 2009). The linear baryon 
power spectrum is related to the DM power spectrum by 
(Fang et al. 1993) 



Pg\k,z) = 



(2) 



where the cut-off length is the scale below which baryon 
density perturbations are smoothed due to physical ef- 
fects like Jeans dissipation or shock heating (Klar & 
Miicket 2010); for the WHIM, shock heating is assumed 
to be the dominating process. At any given redshift the 
comoving scale L cu t is determined by the condition that 
the linear velocity perturbation v(x, z) is equal to or 

larger than the sound speed c s — (kBTiGM(z)/m p ) 1 ^ 2 
at each redshift. Here m p is the proton mass and Tjqm 
is the mean Intergalactic Medium (IGM) temperature. 
At redshifts z < 3, T IGM = 10 3 6 - 10 4 K and its varia- 
tion with redshift is small (Tittley & Meiksin, 2007). In 
addition, L cut ~ Lq(1 + z) 1 / 2 , with Lq a constant and if 
T ~ 10 4 K at present then Lq ~ 1.7/i _1 Mpc. Hereafter, 
Lq will be parameterized by Tjgm (see Suarez- Velasquez 
2013 for details). 



The number density of electrons n e in the IGM can be 
obtained by assuming equilibrium between recombina- 
tion and photo-ionization and collisional ionization. At 
temperatures in the range 10 5 — 10 7 K and density con- 
trasts 6 < 100, the gas can be considered fully ionized. 
The two-point correlation function of the spatial varia- 
tions of the electron pressure is given by 



C(0) 



(Si(x, z)S , 2(i', z'))dzdz'. (3) 



where the integration is along lines of sight separated by 
an angle 0. The TSZ WHIM temperature anisotropy 
is AT = y c G{v), with y c — ko~T /m e c 2 f T e n e dl the 
Comptonization parameter along the line of sight, n e , 
T e the electron density and electron temperature, m e c 2 
the electron annihilation temperature, k the Boltz- 
mann constant, or the Thomson cross section, G(v) = 
(xcoth(a;/2) — 4), with x — hv/kT^ the reduced fre- 
quency, To the CMB temperature, dl the line element 
along the line of sight and dz its corresponding redshift 
interval. Then if S\ = S% = {kaT / m e c 2 )n e T e {dl / dz), 
from eq. ([3]) we can obtain the contribution of the WHIM 
to the power spectrum of CMB temperature anisotropics 
as 

r+l 

C e = 2ir C(0)P t (cos0)dcos0 (4) 



where Pi is the Legendre polynomial of multipole £. In 
Atrio-Barandela & Miicket (2006) we modeled the gas 
as a polytrope. The amplitude of the resulting power 
spectrum was strongly dependent on the polytropic in- 
dex; for some model parameters, it would be larger than 
10Q(fj,K) 2 . Recently, Suarez- Velasquez et al. (2013) have 
refined the model for the shock-heated IGM by deriving 
the relation T e versus n e from a fit to the phase diagrams 
obtained in various hydrodynamical simulations. The fit 
log 10 T e = 8 - 2[log 10 (4 + x 6 )]- 1 with b = a + x' 1 and 
x = n e /n,B the electron density in units of the mean 
baryon density, reproduces well the phase diagram in 
Kang et al. (2005). Then, the model is parameterized by 
the IGM temperature and equation of state with param- 
eters in the range T IGM = [10 3 ' 6 ,10 4 ]K and a = [1,4]. 
For this parameterization, temperature anisotropics grow 
with increasing electron temperature and decreasing cut- 
off length Lq. 

3. CROSS-CORRELATION OF MATTER DENSITY 
TEMPLATES AND CMB MAPS. 

The search of the WHIM contribution to the tem- 
perature anisotropics of the CMB was pioneered by 
Hernandez-Monteagudo et al. (2004). We correlated 
the first year of WMAP lyr with templates of projected 
matter density constructed from the Two Micron All Sky 
Survey (2MASS) galaxy catalog; all significant TSZ con- 
tributions were originated by clusters of galaxies and no 
evidence of the WHIM was found. A second approach 
was tried by Genova-Santos et al. (2009, 2013), who 
used a Monte Carlo Markov Chain to find the contribu- 
tion in the CMB power spectrum but the evidence was 
not statistically significative. In this article we revisit the 
cross-correlation approach using templates of the density 
field reconstructed from 2MRS by Kitaura et al. (2012). 
The reconstruction technique is based on a Bayesian 
Networks Machine Learning algorithm (the KiGEN-code) 



3 



which self-consistently samples the initial density fluc- 
tuations compatible with the observed galaxy distribu- 
tion and a structure formation model given by second 
order Lagrangian perturbation theory (2LPT). We have 
used the Augmented Lagrangian Perturbation Theory 
(ALPT; Kitaura & Hess 2012) to perform constrained 
simulations from the initial conditions found with ki- 
GEN. It improves previous approximations at all scales 
by combining 2LPT with the spherical collapse model 
and shows a higher correlation with the N-body solution 
than previous methods. This approach enables us to find 
non-linear structures like filaments in the cosmic web to 
great accuracy. The number of solutions compatible with 
the observations of a galaxy sample is degenerate due to 
shell crossing and redshift distortions and the method 
provides an ensemble of reconstructions useful to esti- 
mate the uncertainties associated with the technique and 
intrinsic to the data (see Kitaura 2012 for details). Work 
to perform self-consistent reconstructions implementing 
ALPT within the KiGEN-code is in progress (Hess et al. 
2013 in preparation). 

The reconstruction of the matter density field requires 
evaluation of FFTs and, therefore, it is carried out 
on a cubic box. Cubic boxes of side 160/i -1 Mpc and 
180ft. _1 Mpc were used to check that the impact of bound- 
ary effects in the density field is negligible on spheres of 
radius up to 80/i _1 Mpc. To consider even larger volumes 
one would require to model the "Kaiser Rocket effect" in 
the selection function (Branchini et al. 2012) which is be- 
yond the scope of this work. One such reconstruction, de- 
noted by M, is represented in Fig.QJ,. The matter distri- 
bution shows well defined filaments, characteristic for the 
mildly non-linear regime. The sky is represented using 
Healpix (Gorski et al. 2005) with resolution N S id e = 128, 
that corresponds to an angular resolution of 27.5'. We 
used a linear scale saturated at 500 galaxies per pixel 
for better visualization. In Fig. [1}d we represent the dif- 
ference between two reconstructions. While the cosmic 
web displayed in all reconstructions is very similar, the 
extension around very massive structures and the exact 
location of filaments differ from one reconstruction to an- 
other. For example, the uncertainty on the position of 
Coma is 2 — 3/i _1 Mpc, less than 5% of its distance to 
the Local Group. Filaments that are slightly displaced 
appear in the difference map to be running side by side. 
The ensemble of constrained simulations can be used to 
compute the mean and standard deviation of the density 
field in each cell, giving an estimate of the uncertainty 
in the position of the density peaks and filaments. For 
illustration, in Fig[T]; we represent the Wl Differencing 
Assembly of WMAP 7yr data. The galactic and point 
source contaminations are masked using the extended 
temperature analysis WMAP mask KQ75. The data is 
normalized to zero mean and unit variance outside the 
mask. 

In Fig. [5] we represent the mean power spectra of 24 
different template reconstructions of the matter distri- 
bution like in Fig. [T^, that are compatible with 2MRS 
(solid blue line) and the rms dispersion around the mean 
(shaded area). The dot-dashed (green) line represents 
the power spectrum of the projected density of 2MASS 
galaxies used in this reconstruction, corrected by inverse 
weighting with the selection function. To avoid the in- 
stabilities due to small divisors, 37 pixels with number 



density of galaxies equal to or larger than 10 3 galax- 
ies were eliminated. All templates were normalized to 
zero mean and unit variance. Since the projected dis- 
tribution of galaxies is close to a Poisson point process, 
the power spectrum of 2MASS galaxies is roughly con- 
stant and very different, at high multipoles, from the 
spectrum of the reconstructed matter density field. At 
those scales, the matter distribution is smoother and the 
power falls. The solid red line represents the power spec- 
trum of the baryon density distribution computed by 
taking S\ = S% = n e in eq. ([3]), with a cut-off length 
of L ~ 1.13ft. _1 Mpc. We found that the best fit to 
the power spectrum of the reconstructed template was 
obtained when the integration in eq. ((3]) was restricted 
to the interval 60 — 120/i -1 Mpc, while the template of 
Fig. Hk, reproduces the local volume out to 80/i~ 1 Mpc. 
This discrepancy is only apparent. Even if the template 
has been reconstructed up to 80ft. _1 Mpc, it does include 
higher modes since the method uses a 160/i _1 Mpc box. 
Second, we have to take into account sampling variance; 
the local universe lacks power up to 30ft. _1 Mpc (see e.g. 
Courtois et al. 2012). Spikes and oscillations in the spec- 
tra reflect the sample variance associated with observing 
a single universe. While the distribution of the WHIM 
and matter is different from that of galaxies at small 
scales, they coincide at large scales, as it could be expected 
if galaxies traces the overall matter distribution. In the 
inset we represent in a logarithmic scale the power spec- 
trum at low multipoles to show that the reconstructed 
matter density field and the galaxy distribution on those 
scales are very similar. Lines follow the same convention 
that in the main plot. To facilitate the comparison, the 
power spectrum of galaxies in this inset was multiplied 
by a factor of 20 to bring its amplitude closer to that of 
the power spectrum of the reconstructed matter density 
templates. 

Except at very large scales, where the power spec- 
trum is dominated by a few modes and sample variance 
is large, the templates of the reconstructed matter den- 
sity field agree with the log-normal model at all angular 
scales. This agreement demonstrates that templates like 
the one shown in Fig [T^l are very well described by our 
log-normal model, but it does not prove that the WHIM 
is stored in filaments. It only indicates that we can use 
the log-normal model to predict the value of observable 
quantities that can be later compared with the data. To 
this purpose, electron overdensities will be normalized to 
unit variance so the cross-correlation is given in units of 
temperature. For the template of Fig[lji and for model 
parameters in the fiducial ranges Tjom = [10 3,6 , 10 4 ]K 
and a — [1, 4] the correlation at the origin is 0.03— 0.3/^K. 
For a different parameterization that fits the Cen & Os- 
triker (2006) phase diagram, the amplitude would be 
0.01/zK, slightly lower. 

Contributions coming from gas in clusters or in the 
CGM are not included in our formalism. First, the log- 
normal model is restricted to overdensities 5 < 100, while 
in clusters and in the CGM are 5 ~ 500- 1000 (Fukugita 
& Peebles, 2004). Second, the distribution of galaxies or 
clusters is not well described by our log-normal model. 
For instance, the CGM would induce anisotropies with 
the spatial distribution of galaxies and not filaments that, 
as Fig. [2] shows, are very different. 



4 



If a template M is a good tracer of the gas in fila- 
ments, cross-correlation with CMB data will give us an 
estimate of the CMB temperature anisotropics generated 
by the WHIM measured, for example, by WMAP. Since 
the reconstruction is not unique, we need to quantify the 
uncertainty introduced by using a template that does not 
exactly describe the distribution of baryons. To speed up 
the calculation we degraded all our templates to a res- 
olution of ~ 1.8°, corresponding to Healpix N S id e = 32. 
At this resolution, at separation of 1° the correlation 
function is between first neighbors, at 2° between sec- 
ond neighbors, etc. In Fig. [3] we represent the cross- 
correlation of the template M with the 24 different re- 
constructions of the matter density field from 2MRS. 
The solid (black) line represents the autocorrelation func- 
tion of M; the dashed (blue) line represents the mean 
of the cross-correlation of M with the other 23 different 
reconstructions and the shaded area is the rms disper- 
sion about the mean. At the origin, the autocorrelation 
function and the mean of the cross-correlations differ by 
20%. Outside the origin, the autocorrelation is similar to 
the mean and very well within the la error bar, indicat- 
ing that the uncertainty on the position of the density 
peaks in the matter reconstruction of the 2MRS cata- 
log is ~ 2°. Fig. [3] is very illustrative of the accuracy 
of our method. Using resolution N S id e — 32 instead of 
N s ide = 128 we loose information about the correlation 
on scales 0.5° — 2°, but on those scales the reconstruction 
of the density field is uncertain (see Fig [TJd) . If we use a 
reconstruction of the density field that is not fully coin- 
cident with the true distribution of baryons in the local 
universe, we can expect to underestimate the true cor- 
relation function by the 20% but the effect is negligible 
outside the origin. Then, by restricting our analysis to 
Healpix resolution N S id e — 32 not only the computation 
is faster, we also average over the uncertainties on our 
template that only contribute with a 20% uncertainty 
on the correlation at zero lag. With this resolution there 
are no differences between WMAP 7yr and 9yr data and 
the former were used. 

4. APPLICATION TO WMAP DATA 

Since the log-normal model describes reasonably well 
the power spectrum of the matter density then cross- 
correlation of a matter template M with WMAP 7yr data 
would determine the fraction of gas traced by the matter. 
In Fig. |4Jl we represent the cross-correlation of M with 
the eight DAs of WMAP 7yr data. From top to bot- 
tom: solid (blue), dashed (green) and dot-dashed (red) 
lines correspond to the Q, V and W channels. The con- 
tribution of galactic foregrounds and point sources were 
removed using the KQ75 mask that eliminates 27% of 
the sky (see Fig \Vp) . The matter density template M 
was normalized to zero mean and unit variance outside 
the mask. The symmetric solid (black) lines represent 
the rms deviation of the cross-correlations of M with 
1,000 different random realizations of WMAP data that 
include cosmological CMB signal and noise, but no fore- 
ground residuals. The amplitude of the cross-correlation 
at the origin, ~ 1.5/iK, is much larger than the expected 
amplitude of 0.03 — 0.3/^K and it is compatible with be- 
ing due to random alignments between structures in the 
CMB data and in the template. The differences be- 
tween the eight DA are due to the WHIM, noise and 



foreground residuals. Subtracting the correlation of M 
and Wl to the correlation of M with the other seven 
DA removes the contribution due to random alignments 
and yet, due to the frequency variation of the TSZ from 
40 to 90 GHz, it leaves 14-18% of the original TSZ sig- 
nal in the V and Q bands, respectively. If in Fig 
the differences were dominated by TSZ, then the quan- 
tity ((T — Wl)M) = [C vM - C WhM ]/[G(v) - G(Wl)} 
would give the correlation of the template with the map 
of the Comptonization parameter, j/ c -map. In Fig |4j} 
we plot the (M(y c — map))) correlation for WMAP. The 
solid (blue) and dashed (green) lines corresponds to the 
Q and V channels, respectively. The differences between 
the W2, W3, W4 DA with Wl would remove not only 
the intrinsic CMB, but also the foreground and WHIM 
signals. The rms variation of the subtracted correla- 
tions of the W channel provides a simple estimate of 
the error bar due to the instrumental noise, but does 
not account for foreground residuals and sampling vari- 
ance contributions due to the variation of the number of 
pairs with separation angle. The error bars are repre- 
sented in Fig. [4Jd by the symmetrical solid (black) lines. 
The correlations ((Q lj2 - W1)M) and ((V 1)2 - Wl)M) 
are similar and well outside the error bar, but this is not 
evidence of WHIM. In fact, the (M(y c — map)) correla- 
tion should be positive at the origin, not negative. We 
verify that this residual correlation is due to foregrounds 
by masking all CMB data with \b\ < 30°. The results 
are presented in Fig. |4]:. Lines follow the same conven- 
tion as in the panel (b) . By restricting the correlation to 
the polar caps, ((Qi,2 — W1)M) decreased by 30% while 
((^1,2 — W1)M) became positive, with an amplitude of 
0.4/iK at the origin, compatible with noise. Compared 
with Fig. 0J}, the correlation of the Q and V band are 
very different, reflecting the differences in the foreground 
residuals between the Q — Wl and the V — Wl maps and 
an increased sample variance as the fraction of the sky 
removed was larger. 

The previous results indicate that WMAP is not well 
suited to separate the TSZ contribution from that of the 
other components. In this respect, Planck is a much 
better instrument. To forecast its performance, we sim- 
ulated 6 Planck channels with frequencies in the range 
44 — 353GHz. The simulated maps were constrained to 
reproduce WMAP data for multipoles I < 256. The 
maps were downgraded to resolution N S id e — 32 so dif- 
ferences in beamwidth and beam asymmetries have an 
unmeasurable effect. The simulated data contain CMB 
and homogeneous white noise but did not contain fore- 
grounds, foreground residuals or (1/f) noise. At each 
frequency we added a WHIM component following the 
matter distribution in M assuming T e oc n e with a mean 
Comptonization parameter (y c — map) = 0.1//K. In the 
frequency range considered, the amplitude of TSZ effect 
varies from —1.9 at 44GHz to 2.2 at 353 GHz and passes 
through the TSZ null at 217GHz. This channel was later- 
used to subtract the correlation of M with the CMB data 
due to random alignments. In Fig.|4ji we plot the results 
(solid black line), that are very close to each other and 
close to the theoretical expectation, represented by the 
dashed (blue) line. The differences between the theoreti- 
cal curve and the simulated data are due to masking and 
pixelization. If foreground and noise inhomogeneities are 
not important and the template traces the gas distribu- 
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tion, Planck will measure the WHIM TSZ correlation 
very accurately. 

Together with the uncertainties in the reconstruction of 
the density field, there are additional uncertainties asso- 
ciated with how well the galaxy template traces the elec- 
tron pressure. The WHIM could be made of clumps, with 
typical sizes 100/i _1 Kpc like the CGM. At 60/i _1 Mpc, a 
clump of this size would subtend 6'. The 217 and 343GHz 
Planck channels have the largest resolution of the instru- 
ment, 5'; then, in Planck data we can neglect the effect 
of gas inhomogeneities at scales below 100ft. _1 Kpc. Also, 
in our analysis we have assumed a smooth distribution 
of baryons and temperatures within the filaments. As 
an average of the different Kang et al. (2005) models we 
assumed the temperature scaling as T e oc n e so the elec- 
tron pressure in the template scaled as n e T e oc n\. The 
simulations of individual filaments carried out in Klar 
& Mucket (2010, 2012) have characterized the state of 
the gas and its evolution inside the filament in terms 
of the length of the initial perturbation. Their simu- 
lations indicated the existence of multiple phases and 
temperatures in the WHIM. To estimate the uncertainty 
associated with an inhomogeneous distribution of gas 
and temperature within filaments, we repeat our sim- 
ulations but added to the CMB data a TSZ contribu- 
tion with two different temperature scalings in the range 
S G [1,100]: (1) an isothermal gas T ~ const and (2) a 
clumpier temperature distribution, T e oc n\. In both 
cases, the mean amplitude of the y c -map was 0.1/iK. 
These temperature-density relations can be considered 
extreme cases of the temperature variation in the IGM. 
Next, we cross correlate the CMB data with our template 
M, where the electron pressure scales like n e T e oc n\. 
The resulting correlations are represented by dot-dashed 
(red and violet) lines in Fig0Jl. The comparison of these 
cross-correlations with our previous result shows that, at 
N S ide — 32 resolution, the shapes are similar and only 
differ by ~5-10% in the amplitude at the origin. In sum- 
mary, the uncertainty on the cross-correlation of tem- 
plate and CMB data due to (1) how the template traces 
the gas distribution and (2) how the electron pressure 
scales with density within the template add at most an 
uncertainty of 30% at the origin, being almost negligible 
on scales above 2°. 

The final error bar will certainly have non-zero con- 
tributions from foreground residuals, whose distribution 
and amplitude change with frequency, and instrumental 
(1//) noise that will complicate the detection. However, 
first and foremost, the TSZ has a distinctive signature, 
different from that of any other foreground. It changes 
sign from the Rayleigh- Jeans to the Wien part of the 
spectrum. This change of sign will also be present in 
the correlation function, and could be sufficient to de- 
tect the effect. Second, foregrounds residuals are proba- 
bly largest close to the galactic plane, so computing the 
cross-correlation at different galactic latitudes could help 
to disentangle the different contributions, as for WMAP 
data. Third, in our Planck simulations, the template 
contributes with just 0.1/iK. According to our model, 
this is a reasonable expectation for a large fraction of 
the parameter space. Naturally, the contribution due a 
larger volume will be higher so templates constructed 
from deeper galaxy surveys would give statistically more 



significant detections. 

5. CONCLUSIONS. 

We have computed the amplitude of the correlation be- 
tween the WHIM with CMB data. We have shown that 
templates of projected matter density that describe the 
weakly non-linear regime are well represented by our log- 
normal model. This permits us to predict that the am- 
plitude of the cross-correlation of the WHIM with CMB 
data is 0.03— 0.3uK depending on model parameters. As- 
suming that the reconstructed matter density templates 
trace the electron distribution, we computed the cross- 
correlation of the matter density field within a volume of 
160/i _1 Mpc reconstructed using the 2MRS galaxy sur- 
vey. The cross-correlation with WMAP data showed it to 
be dominated by random alignments of CMB structures 
and galaxy filaments. The differences between DA could 
be due to the WHIM signal, but also to noise and fore- 
ground residuals. We checked that foreground residuals 
are the most likely source of the correlation by showing 
the results varied significantly when removing data with 
|6| < 30°. 

Using 2MASS galaxies and WMAP data, Lavaux et al. 
(2013) described evidence that baryons are distributed 
in galactic coronae. This result is in apparent contra- 
diction with Hernandez-Montegudo et al. (2004) where 
cross-correlation of WMAP lyr data with templates of 
projected galaxy density showed no evidence of ionized 
gas outside known clusters of galaxies and, to some de- 
gree, with the results presented here. Even if our matter 
density templates are different from a pure galaxy tem- 
plate their power spectra overlap (see Fig [2]) and any sig- 
nal traced by galaxies must also be traced at some level 
by our templates. In addition, the contribution due to 
random alignments, that we have shown to be dominant 
for WMAP, was not quantified. Due to the difference in 
methodology, our results are not fully comparable and 
consistency between both methods will increase the sta- 
tistical significance of any detection. 

Planck, with its large frequency coverage and high res- 
olution is a more adequate instrument than WMAP to 
search for the WHIM contribution. First, the 217GHz 
channel is very close to the TSZ null, allowing to remove 
spurious correlations that dominates the signal. Second, 
the instrument contains measurements on the Rayleigh- 
Jeans and Wien part of the spectrum. If the correlation is 
due to the WHIM, it will vary with frequency and change 
sign. Reconstructions of the matter density field using 
deeper surveys like the Sloan Digital Sky Survey will al- 
low to probe the baryon distribution to higher depths, 
enhancing the signal. Alternatively, the reconstructed 
density field allows to select filaments aligned with the 
line of sight, where the optical depth of the WHIM would 
be larger. The higher resolution of Planck could facilitate 
to identify the WHIM signal in those regions and clar- 
ify its spatial distribution, whether it is distributed as a 
network of filaments or is stored in the galactic coronae. 
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(a) ONE 2 MASS DENSITY FIELD RECONSTRUCTION 
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(b) DIFFERENCE OF TWO RECONSTRUCTIONS 
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Fig. 1. — Templates and CMB data at Healpix resolution N ai d e = 128. (a) Template M of projected matter density re- 
constructed from the 2MASS catalog, showing the filamentary structure characteristic of the weakly non-linear regime, (b) 
Difference of two reconstructed density fields showing small displacements in the location of filaments and on the extent of 
massive structures, (c) WMAP data of the Wl Differencing Assembly normalized to zero mean and unit variance; the galactic 
plane has been masked using the KQ75 galactic and point source mask. 
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Fig. 2. — Mean power spectra of the 24 reconstructed matter distribution templates including the one represented in FigQJ, 
(solid blue line) and rms dispersion around the mean (shaded area). The dashed (green) line is the power spectrum of the density 
of 2MASS projected galaxies corrected by the selection function and the solid (red) line is the spectrum of baryon distribution 
given by our log-normal model. All power spectra were normalized to zero mean and unit variance. In the inset we represent 
the low order multipoles in a log-scale; the amplitude of the power spectrum of galaxies has been multiplied by a factor 20 to 
facilitate the comparison. 
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Fig. 3. — Cross-correlation of the different reconstructions of the matter density field. The autocorrelation of the template 
M of Fig is represented by the solid (black) line while the blue dashed line and shaded are represents the mean and rms 
dispersion of the cross-correlation of M and the other 23 different reconstructions. For convenience, the values on the X-axis 
have been shifted by one degree. 
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Fig. 4. — (a) Cross correlation of the template M, represented in Fig. with the 8 DA of WMAP 7yr data. From top to 
bottom, the solid (blue), dashed (green) and dot-dashed (red) lines correspond to the Q, V and W channels. The symmetric solid 
black lines correspond to the rms cross-correlation of the template with 1000 randomly generated Wl maps, (b) Cross-correlation 
of the Q (solid, blue) and V (dashed, green) channels after subtracting the Wl cross-correlation, renormalized to account for the 
removed TSZ contribution, i.e., {(T-W1)M) = [C„(0)-Cwi(0)]/[G(v)-G(Wl)]. Black lines represent the rms of the W2, W3, 
W4 DA after subtracting the Wl cross-correlation, (c) As in (b) but the CMB-template correlation was computed only on pixels 
with \b\ > 30°. (d) Planck forecast: Solid (black) lines represent ((T -T 2 n G Hz)M) = [C„(0) - G-xyighz (0)]/[G(v) -G{2Y!GHz)] 
for the 44-353GHz channels; the dashed (blue) line is the log-normal prediction. The dot-dashed lines (red and violet) represent 
the cross-correlations of the template M with Planck simulated data containing a TSZ WHIM contribution with T e oc const, n\ 
in the range 8 S [1, 100], respectively. Like in Fig[3]the x-axis has been shifted by one degree. 
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